AD-A107  864  MASSACHUSETTS  INST  OF  TECH  CAMBRID6E  DEPT  OF  MATERIA— ETC  F/G  20/2 

ELECTRONIC  STRUCTURE  OF  GRAIN  BOUNDARIES  AND  INTERFACES  IN  POLY~ETC<U> 
NOV  81  M  H  SUKKAR*  H  L  TULLER •  K  H  JOHNSON  N00014-81-K-0499 
unclassified  TR-2  nl 


FILE  copy  AD  A1 078 6 4 


^ECUrtlTV  CLASSIFICATION  OF  THIS  PAGE  fWhmn  Dm •  En(«r«dN 


REPORT  DOCUMENTATION  PAGE 


.  REPORT  NUMBER 


.uia 


*.  TITLE  (md  Submit) 

Electronic  Structure  of  Grain  Boundaries  and 
Interfaces  in  Polycrystal  line  Zinc  Oxide 


READ  tSSTRUCTlONS 
BEFORE  COMPLETING  FORM 


S.  TYPE  OF  REPORT  «  PERIOD  COVERED 


Interim 


ft.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AyTHORf*J 

M.  H.  Sukkar,  H.  L.  fuller,  and  K.  H.  Johnson 


*■  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


I.  CONTRACT  or  grant  NUMBERCR) 


N00014-81 -K-G499 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  *  WORK  UNIT  NUMBERS 


Task  No.  NR  056-757 


12.  REPORT  DATE 

November  18,  1981 


<<.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Office  of  Naval  Research 
Department  of  the  Navy 

Arlington,  Virginia  22217  _ 


.  MONITORING  AGENCY  name  ft  ADDRESSflf  dlilerent  from  Controlling  OtUca)  IS.  SECURITY  CLASS,  (ot  I  e  report) 


13.  NUMBER  OF  PAGES 


13a.  DECLASSI  FI  CATION /DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (of  thla  Raport) 


Approval  for  public  release;  distribution  unlimited. 


19.  KEY  WORDS  (Contlnum  on  reeorae  a  I  da  it  nacaaaary  and  Identity  by  block  number) 


Electronic  structures;  grain  boundaries;  polycrystall ine;  zinc  oxide. 


20.  ABSTRACT  (Contlnua  on  reverae  aide  It  necaaaory  mnd  Idantity  by  block  number) 

'Preliminary  theoretical  models  for  the  electronic  structure  of  grain 
boundaries  and  interfaces  in  polycrystalline  ZnO  have  been  constructed  on  the 
basis  of  self-consistent-field  X-alpha  scattered-wave  (SCF-Xot-SW)  cluster 
molecular-orbital  calculations.  The  disposition  and  character  -if  the  inter¬ 
face  states,  relative  to  the  valence  and  conduction  bands  of  the  otherwise  per¬ 
fect  crystalline  material,  have  been  studied  for  clusters  representing 
coordinatively  unsaturated  Zn  surface  sites  and  molecular  0 2  chemisorption 

DO  ,  j2U  7J  1473  EDITION  OF  I  NOV  •»  It  OMOLtTC  rjn  1  f-jgH 

S/N  01 02- LF -01 4-6601  -  -  _ 


Unclassified 

SECURITY  CLASSIFICATION  OF  TMIS  PAGE  fWliwn  Data  «n land) 

>  8  1  11  24  166 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whrni  Dim  Kntmd) 


thereon.  The  possible  effects  of  the  resulting  interface  states  on  electron 
transport  at  grain  boundaries  in  ZnO  varistors  have  been  addressed. 


UNCLASSIFIED  ,  _ _ 

SECURITY  CLASSIFICATION  OF  TI»S  FAOEf»*i»  Omm  SntarMQ 


OFFICE  OF  NAVAL  RESEARCH 


Contract  N00014-81-K-0499 
Task  No.  NR  056-757 


TECHNICAL  REPORT  NO.  2 


Electronic  Structure  of  Grain  Boundaries  and  Interfaces  in  Polycrystalline 
Zinc  Oxide 


by 


M.  H.  Sukkar,  H.  L.  Tuller,  and  K.  H.  Johnson 
Department  of  Materials  Science  and  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 


November  18,  1981 

Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 
Approved  for  Public  Release:  Distribution  Unlimited 
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ABSTRACT 

Preliminary  theoretical  models  for  the  electronic  structure 
of  grain  boundaries  and  interfaces  in  polycrystalline  ZnO 
have  been  constructed  on  the  basis  of  self-consistent- field 
X-alpha  scattered-wave  (SCF-Xo-SW)  cluster  molecular-orbital 
calculations.  The  disposition  and  character  of  the  inter¬ 
face  states,  relative  to  the  valence  and  conduction  bands 
of  the  otherwise  perfect  crystalline  material,  have  been 
studied  for  clusters  representing  coordinatively  un¬ 
saturated  Zn  surface  sites  and  molecular  02  chemisorption 
thereon.  The  possible  effects  of  the  resulting  inter¬ 
face  states  on  electron  transport  at  grain  boundaries  in 
ZnO  varistors  have  been  addressed. 


INTRODUCTION 

Electronic  ceramics  whose  properties  are  largely  controlled  by  grain 
boundary  phenomena  have  received  a  good  deal  of  attention  in  recent  years.  One 
such  device,  the  ZnO-based  varistor,  which  is  resistive  at  low  voltages,  "breaks 
down"  at  higher  voltages  (hundreds  of  volts/cm)  to  give  a  highly  nonehmic 
regime  in  which  the  current-voltage  relation  is  described  by  a  power  law, 

I  =  kV“  (a  >  1).  Early  models  [1J  attempted  to  ascribe  the  nonlinear  varistor 
action  to  the  existence  of  second-phase  layers  of  insulating  material  separat¬ 
ing  the  semiconducting  grains.  Recent  high  resolution  transmission  electron 
microscopy  [2]  shows,  however,  this  not  to  be  the  case.  Although  a  considerable 
number  of  alternative  models  have  been  proposed  in  recent  years  [3,4]  they  all 
propose  in  common  the  existence  of  a  charged  grain  boundary  interface  bounded 
by  space  charge  layers  penetrating  into  the  adjoining  grains.  All  models 
either  implicitly  or  explicitly  rely  on  the  existence  of  low-lying  interface 
states  or  traps  at  the  intergranular  layers  to  provide  for  the  formation  of  the 
potential  barrier  at  the  grain/grain-boundary  interface.  In  particular,  the 
recent  models  of  Pike  and  Seager  [5]  quantitatively  derive  expressions  for  the 
temperature  and  voltage  dependence  of  the  potential  barrier  with  the  density- 
of-states  function  of  the  Interfacial  traps  serving  as  a  key  parameter.  The 
purpose  of  our  investigation  is  to  examine  the  nature  of  the  electronic  struc¬ 
ture  of  ZnO  surfaces  and  Interfaces  from  first  principles,  and  subsequently  to 
attempt  to  relate  our  findings  to  experimental  observations.  This  paper  serves 
as  our  first  effort  in  this  direction. 

Zinc  oxide  which  forms  in  the  wurtzite  structure  possesses  a  band  gap  energy 
of  3.3eV  at  room  temperature.  It  can  nevertheless  be  made  highly  semiconduct¬ 
ing  (n-type)  either  by  heating  in  a  reducing  atmosphere  to  obtain  an  excess  of 
Zn  or  by  doping  with  donor  Impurities  (e.g.  At).  Under  these  conditions  donor 
ionization  energies  are  reported  to  be  low  (,v<0.05eV),  while  electron  mobilities 
are  typically  ^  200  cnr/V-sec  at  300°K.  Acceptor  impurities  (e.g.  Li),  on  the 
other  hand,  can  drop  the  conductivity  of  these  materials  many  orders  of  magni¬ 
tude.  Likewise,  annealing  ZnO  in  oxygen  can  be  expected  to  have  a  similar 
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effect  [6].  In  the  following  analysis  we  examine  in  some  detail  the  effects  of 
excess  oxygen  on  ZnO  interfaces  but  limit  our  discussion  to  the  nature  of  the 
interface  states  and  how  they  are  modified  upon  chemisorption  of  oxygen. 

THE  SCF-Xa-SW  CLUSTER  MOLECULAR-ORBITAL  METHOD 

The  inherently  local  nature  of  the  ZnO  grain  boundary  problem,  coupled  with 
the  importance  of  chemical  bonding  (e.g.  oxygen  chemisorption)  at  the  grain 
boundary  interface,  make  traditional  reciprocal  (fc)-space  band-theoretical 
methods  of  solid  state  physics  both  Inappropriate  and  impractical  for  determin¬ 
ing  grain-boundary  electronic  structure/properties  relationships.  Real-space 
molecular-orbital  theory  is  more  appropriate  for  this  problem  because  it 
focuses  on  the  electronic  structure  of  a  small  cluster  of  atoms  which  simulates 
the  local  environment  of  the  grain  boundary.  The  feasibility  of  this  approach 
is  shown  by  the  successful  modeling  of  sulfur  embrittlement  of  nickel  grain 
boundaries  by  Messmer  et  al  [7]  using  the  SCF-Xa-SW  cluster  molecular-orbital 
method  [8J. 

The  theoretical  formalism  and  computational  procedure  underlying  the  SCF-Xa- 
SW  molecular-orbital  method  have  been  thoroughly  documented  elsewhere  [8].  It 
is  sufficient  to  summarize  here  that  in  applying  this  technique  to  grain  boun¬ 
daries  in  ZnO,  appropriate  clusters  of  atoms  are  chosen  which  are  models  for  the 
local  arrangement  of  atoms  at  the  grain  boundary  interface.  Each  cluster  is 
partitioned  into:  I.  atomic,  II.  interatomic,  and  III.  extramolecular  regions  of 
spherically  and  volume  averaged  potentials.  Including  the  Xa  approximation  to 
electron-electron  exchange  and  correlation,  based  initially  on  free-atom  and 
free-ion  Xa  potentials.  Schrbdinger's  equation  is  then  solved  numerically  in  a 
composite  partial-wave  representation  throughout  regions  I,  II,  and  III,  and  the 
solutions  along  with  their  first  derivatives  are  joined  continuously  via  scat¬ 
tered-wave  theory  across  the  spherical  boundaries  separating  these  regions. 

These  continuity  conditions  determine  the  cluster  molecular-orbital  wavefunc- 
tions,  and  their  associated  energy  levels,  which  are  then  occupied  with 
electrons  according  to  the  Pauli  principle  and  Fermi  statistics.  From  the 
occupied  orbital  wavef uncti ons ,  a  new  electronic  charge  density  is  calculated 
which,  in  turn.  Is  used  to  construct  a  new  potential  throughout  the  cluster, 
and  Schrodinger's  equation  is  again  solved  for  a  new  set  of  molecular-orbital 
wavef uncti ons  and  energy  levels.  This  entire  process  is  iterated  until  self- 
consistency  of  the  electronic  energy  levels  and  charge  distributions  is  attained. 
The  procedure  has  been  programmed  for  rapid  computation  on  a  digital  computer. 

PRELIMINARY  THEORETICAL  MODEL 

The  rectifying  behavior  of  ZnO  surfaces  (electroded  with  Pt)  exposed  to 
oxygen,  which  disappears  when  heated  to  ^600°K  in  a  vacuum  and  reappears  when 
O2  is  reintroduced  [9],  strongly  suggests  a  direct  correlation  between  excess 
oxygen  at  ZnO  interfaces  and  nonlinear  I-V  characteristics.  This  correlation 
Is  further  supported  by  the  disappearance  of  the  nonlinear  I-V  characteristics 
of  ZnO  varistor  materials  upon  annealing  in  a  reducing  atmosphere  [10]. 

In  order  to  provide  a  theoretical  basis  for  this  correlation,  preliminary 
SCF-Xa-SW  molecular-orbital  calculations  of  the  type  described  in  the  precee- 
dlng  section. have  been  carried  out  respectively  for  molecular  02  and  Zn03  4- 
and  0?Zn03  clusters.  The  O2  molecule  calculations  were  carried  out  in  order 
to  (1)  estimate  the  positions  of  chemisorbed  02  levels  with  respect  to  the  bulk 
ZnO  levels  and  (2)  enable  a  comparison  of  the  occupation  of  its  energy  levels 
before  and  after  chemisorption  on  ZnO  surfaces.  The  ZnOj  4-cluster  was  chosen 
as  an  initial  model  for  ZnO  free  surfaces  because  of  its  simplicty  ana  its 
preservation  of  underlying  bulk  character.  This  cluster  may  be  thought  of  as 
representing  either  (1)  a  coordinatively  unsaturated  Zn  site  on  a  (0001)  Zn 


Fig.  2.  Contour  map  for  the  8e, 
molecular-orbital  wavefunctlon  Af 
a  ZnO^"  cluster.  The  solid  and 
dashed  contours  represent  opposite 
phases  of  the  wavefunctlon. 


Fig.  3.  Contour  map  for  the  9e, 
molecular-orbital  wavefunction  Af  an 
09Zn0-4"  cluster.  The  solid  and  dash¬ 
ed  contours  represent  opposite  phases 
of  the  wave function. 


surface,  (2)  a  step  site  on  a  {10T0}  prismatic  surface,  (3)  an. oxygen  vacancy  on 
a  (OQOT)  0  surface,  or  (4)  a  bulk  oxygen  vancancy.  The  O-ZnO.*”  cluster  repre¬ 
sents  an  02  molecule  chemisorbed  on  the  coordinatively  unsaturated  Zn03 
cluster,  with  the  02  molecular  axis  perpendicular  to  the  plane  formed  By  the 
three  oxygens.  This  choice  of  02. orientation  was  made  to  preserve  the  C3 
point  group  symmetry  of  the  Zn03  "  cluster,  although  preliminary  calculations 
for  02  oriented  "parallel"  to  the  surface  indicate  that  the  main  conclusions 
about^the  effects  of  02  chemisorption  on  interface  electronic  structure  are  un¬ 
altered.  These  clusters  are  not  large  enough  to  distinguish  between  various 
ZnO  surfaces  or  between  unsaturated  Zn  surface  sites  on  a  (0001)  surface  and 
surface  oxygen  vacancies  or  steps.  Nevertheless,  it  is  felt  that  the  study  of 
the  electronic  structures  of  these  simpler  clusters  is  a  necessary  first  step 
toward  understanding  grain  boundary  electronic  structure  and  oxygen  incorpora¬ 
tion  therein.  The  resulting  cluster  energy  levels  are  shown  in  Fig.  1,  where 
they  are  compared  directly  with  the  SCF-Xo  orbital  energies  for  the  02  mole-' 
cule  and  Zn  and  0  atoms.  The  levels  are  labeled  according  to  the  irreducible 
representations  of  the  appropriate  (e.g.  C3  )  synmetry  groups,  and  the  highest 
occupied  levels  (the  "Fermi  levels")  at  absolute  zero  temperature  are  indicated 
by  horizontal  arrows.  The  electron  spin  occupancy  of  the  CL  molecular-orbital 
energy  levels  is  indicated.  Also  shown  schematically  for  comparison  is  the 
bulk  ZnO  band  structure,  which  has  been  successful ly  modeled  by  previous 
SCF-Xa-SW  molecular-orbital  calculations  for  a  ZnO.0”  cluster  representing  the 
nearest-neighbor  environment  in  the  bulk  material  [11].  The  coincidence  of  the 
Fermi  level  with  the  bottom  of  the  conduction  band,  as  indicated  by  the  arrow, 
denotes  the  usual  n-type  material. 

It  is  found  that  a  coordinatively  unsaturated  Zn  surface  site  (or  surface  oxy¬ 
gen  vacancy),  as  simulated  by  the  Zn03  cluster,  produces  occupied  states  in 
the  lower  part  of  the  otherwise  unoccupied  Zn  s,p-like  conduction  band  of  in¬ 
trinsic  bulk  ZnO.  The  predominantly  Zn  4p-like  character  of  the  highest  oc¬ 
cupied  of  these  conduction-band  states,  spatially  directed  away  from  the  sur¬ 
face  Zn  site  and  antibonding  with  respect  to  the  02p  orbitals,  is  clearly  vis¬ 
ible  in  the  wavefunction  contour  map  for  the  thermally  occupied  (T~  300°K)  8e, 
orbital  of  the  Zn03  cluster,  plotted  in  the  X-Z  plane  of  the  Zn  atom  and  a 
nearest-neighbor  O^atom.  Because  of  the  embedding  of  these  surface  states  in 
the  conduction  band,  they  function  as  electron  donors,  consistent  with  the 
n-type  behavior  of  oxygen  vacancies  in  ZnO. 

Along  with  the  s,p-like  surface  states  in  the  conduction  band,  the  anti- 
bonding  Zn(3d)-0(2p)  character  in  the  otherwise  primarily  0(2p)-like  valence 
band  is  significantly  enhanced  by  the  coordinative  unsaturation.  This  en¬ 
hanced  Zn(3d)  character  in  the  valence  band  can  also  be  interpreted  as  surface 
states  and  may  be  a  major  contribution  to  the  well  known  catalytic  activity  of 
ZnO  surfaces.  The  presence  of  surface  states  in  both  the  conduction  and  valence 
band  and  absence  of  such  states  in  the  band  gap  are  consistent  with  previous  ex¬ 
perimental  [12-14]  and  theoretical  [15]  findings. 

The  calculated  electronic  structures  of  QJ  before  and  after  chemisorption  on 
a  ZnO  surface,  as  represented  by  the  molecular-orbital  energies  of  02  and 
02Zn03  in  Fig.  1(a)  and  1(b),  respectively,  indicate  that  there  is'Mndeed  a 
net  transfer  of  approximately  one  electron  to  the  02  molecule  upon  chemisorp¬ 
tion,  as  suggested  In  the  literature  [16].  Over  the  temperature  range, 

300°K  ^T  %  700°K,  for  which  non-dissociative  02  chemisorption  (as  distinguished 
from  02  physisorption  for  T  *  300°K)  on  ZnO  is  Believed  to  occur  [16],  the  pre¬ 
sent  theoretical  model  suggests  that  thermally  populated  "surface  states"  em¬ 
bedded  in  the  conduction  band  [e.g.  the  Qe-\  orbital  in  Fig.  1(c)  for  the  re¬ 
presentative  Zn03  "  cluster],  and  having  the  spatially  localized  Zn  4pir*  cha¬ 
racter  shown  in  Fig.  2 .effectively  transfer  one  electron  per  coordinatively  un¬ 
saturated  Zn  surface  site  to  each  02  molecule  through  symmetry  conserving  over¬ 
lap  with  the  half  empty  02  lir  molecular  orbital  shown  in  Fig. 1(a).  This  elec- 
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tron  transfer,  followed  by  considerable  orbital  rehybridization  within  the 
valence  band  at  the  O^-ZnO  interface,  results  In  an  empty  conduction  band  and 
partially  occupied  02  **-11ke  states  coinciding  with  the  top  of  the  valence 
band.  .The  latter  states  are  represented  by  the  9e, , orbital  of  the  model 
O^nO,*-  cluster  In  Fig.  1(b)  having  the  orbital  topology  shown  In  Fig.  3. 

Note  also  In  Fig.  1(b)  the  one-to-one  correspondence  between  the  various  02 
bonding  molecular-orbital  energy  levels  below  the  lir  level  and  energy  levels, 
of  the  02Zn0,^"  cluster  lying  well  below  the  top  of  the  valence  band,  consis¬ 
tent  witn  the  assumption  of  non-dissoclatlve  02  chemisorption  on  ZnO.  The 
presence  of  partially  empty  02  **-like  states  6f  the  type  shown  in  Figs.  1(b) 
and  3  at  the  top  of  the  valence  band  implies  p-type  ("hole")  character  at  the 
0?-Zn0  interface,  in  distinct  contrast  to  the  n-type  ("electron")  character  of 
the  bulk  and  surface  electronic  structures  of  ZnO  shown  in  Figs.  1(d)  and  1(c), 
respectively.  This  dramatic  change  from  n  to  p  character  when  occupied  elec¬ 
tronic  states  in  the  conduction  band  are  compensated  for  by  02  chemisorption, 
assuming  continuity  of  the  Fermi  energy,  reflects  the  Schottky-barrier-like 
"band  bending"  that  is  believed  to  occur. at  the  grain  boundaries  in  ZnO  var¬ 
istors  [3-5].  The  possible  role  of  the  02  **-like  hole  states  at  the  top  of  the 
valence  band  at  an  0?-Zn0  interface  as  intergranular  electron  "trapping"  states 
responsible  for  the  rapid  increase  of  electric  current  with  applied  voltage  in 
ZnO  varistors  [3-5]  is  currently  under  theoretical  investigation  using  the 
SCF-Xo-SW  cluster  molecular-orbital  method. 
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